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plus Structured Interference-Part I: 

On the Maximal Invariant Statistic 
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Abstract —This paper deals with the problem of adaptive 
multidimensional/multichannel signal detection in homogeneous 
Gaussian disturbance with unknown covariance matrix and 
structured deterministic interference. The aforementioned prob¬ 
lem corresponds to a generalization of the well-known Gener¬ 
alized Multivariate Analysis of Variance (GMANOVA). In this 
first part of the work, we formulate the considered problem 
in canonical form and, after identifying a desirable group of 
transformations for the considered hypothesis testing, we derive 
a Maximal Invariant Statistic (MIS) for the problem at hand. 
Furthermore, we provide the MIS distribution in the form of 
a stochastic representation. Finally, strong connections to the 
MIS obtained in the open literature in simpler scenarios are 
underlined. 

Index Terms —Adaptive Radar detection, CFAR, Statisti¬ 
cal Invariance, Maximal Invariants, Double-subspace model, 
GMANOVA, coherent interference. 


I. Introduction 
A. Motivation and Related Works 

T HE PROBLEM of adaptive detection of targets embedded 
in Gaussian interference is an active research field which 
has been object of great interest in the last decades. Many 
works appeared in the open literature, dealing with the design 
and performance analysis of several detectors handling many 
specific detection problems (the interested reader is referred 
to m and references therein). 

It can be shown that most of the aforementioned models can 
be seen as special cases of the model considered by Kelly and 
Forsythe 0, which is very general and encompasses point¬ 
like and extended targets as special instances. The considered 
model allows for training samples which contain random 
interference modeled as an unknown covariance matrix that 
accounts for both clutter and thermal noise, with the implict 
assumption that signal plus noise and noise-only vector sam¬ 
ples share the same covariance matrix, thus determining a so- 
called homogeneous environment. 

The signal model considered in the aforementioned report is 
the well-known Generalized Multivariate Analysis of Variance 
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(GMANOVA) in statistics literature 0, also referred to as a 
“double-subspace” signal model (see for example 01, 0). The 
standard GMANOVA model was first formulated by Potthoff 
and Roy 0 and consists in a generic patterned mean problem 
with a data matrix whose columns are normal random vectors 
with a common unknown covariance matrix. The GMANOVA 
model was later studied in more detail in Q, where maximum 
likelihood estimates of unknown parameters were obtained. 
For a detailed introduction to estimation and detection in 
GMANOVA model (along with few interesting application 
examples) the interested reader may refer to the excellent 
tutorial 0. 

Differently, in this paper we will study a modified version 
of GMANOVA with respect to its classical formulation 0, 
referred to as I-GMANOVA in what follows. The considered 
model allows for the presence of a structured (partially known) 
non-zero mean under both hypotheses. Such disturbance is 
collectively represented as an unknown deterministic matrix, 
which determines an additional set of nuisance parameters 
for the considered hypothesis testing (i.e., other than the 
covariance matrix). The aforementioned model easily accounts 
for the presence of structured subspace interference affecting 
the target detection task. Thus it is clear that taking such inter¬ 
ference into consideration enables the application of this model 
to adaptive radar detection; for instance it may accomodate the 
presence of multiple pulsed coherent jammers impinging on 
the radar antenna from some directions. 

Although several different decision criteria can be consid¬ 
ered to attack composite hypothesis testing problems 0, Go), 
an elegant and systematic way consists in resorting to the so- 
called Principle of Invariance CD, 110. Indeed, the aforemen¬ 
tioned principle, when exploited at the design stage, allows 
to focus on decision rules enjoying some desirable practical 
features. The preliminary step consists in individuating a suit¬ 
able group of trasformations which leaves the formal structure 
of the hypothesis testing problem unaltered. With reference to 
radar adaptive detection, the mentioned principle represents an 
effective tool for obtaining a statistic which is invariant with 
respect to the set of nuisance parameters, therefore constituting 
the basis for Constant-False Alarm Rate (CFAR) rules. Indeed, 
every invariant decision rule can be written in terms of the 
maximal invariant statistic. Therefore, with reference to I- 
GMANOVA model, the principle of invariance allows for 
imposing CFARness property with respect to the clutter plus 
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noise (disturbance) covariance matrix and the jammer location 
parameters. 

It is worth remarking that the use of the invariance principle 
for generic composite hypothesis testing problems IfTTI . IflOl 
(and, more specifically, in the context of radar adaptive detec¬ 
tion) is not new. Indeed, starting from the seminal paper im, 
many works focused on adaptive radar detection problems 
with the use of invariance theory. For example, in ED, 131, 
El, invariance theory was exploited to study the problem 
of single-subspace (adaptive) detection of point-like targets. 
Later, similar works appeared in the open literature dealing 
with the case of a target spread among more range cells 
m. ezi. More recently, the same statistical tool has been 
employed to address the problem of adaptive single-subspace 
detection problem (of point-like targets) in the joint presence 
of random and subspace structured interference in 1181 . In this 
respect, we build upon the aforementioned results in order to 
develop an exhaustive study for the considered I-GMANOVA 
model under the point of view of the invariance. 


B. Summary of the contributions and Paper Organization 

The main contributions of the first part of the present study 
are summarized as follows: 


• We first show that the problem at hand admits a more 
intuitive representation, by exploiting a canonical form 
representation. Such representation helps obtaining the 
maximal invariant statistics and gaining insights for the 
problem under investigation; 

• The group of transformations which leaves the problem 
invariant is identified, thus allowing the search for a MIS. 

• Given the aforementioned group of transformations, the 
canonical form is exploited in order to obtain the MIS, 
which, for the I-GMANOVA model is represented by 
two matrices which compress the original data. Such 
result can be interpreted as the generalization of the 
two-components scalar MIS obtained in the classical 
references na. ca. 

• A theoretical performance analysis of the MIS is ob¬ 
tained, in terms of its distribution. Even though in the 
considered setup the MIS does not generally admit an 
explicit expression for its probability density function 
(pdf), a simpler form of the statistic distribution, by 
means of a suitable stochastic representation, is provided. 

• Finally, the obtained MIS expression is compared with 
similar findings obtained in the literature for simpler 
scenarios, thus showing that the aforementioned cases can 
be seen as special instances of the obtained MIS. 


The explicit expression of the MIS obtained in this first part 
is then exploited to show CFARness of all the detectors 
considered in part II of this work. 

The remainder of the paper is organized as follows: in 
Sec. [U] we introduce the hypothesis testing problem under 
investigation; in Sec. Ill we describe the desirable invariance 
properties and derive the MIS. Sec. [TV] is devoted to the 
statistical characterization of the MIS, while in Sec. [V] we 
particularize the MIS to specific instances and compare it 
with previously obtained results in the open literature. Some 


concluding remarks and future research directions are given 
in Sec. |VI[ finally, proofs and derivations are confined to the 
Appendices. 

Notation - Lower-case (resp. Upper-case) bold letters denote 
vectors (resp. matrices), with a n (resp. A„ jm ) representing 
the nth (resp. the (n, m)th) element of the vector a (resp. 
matrix A); R w , C N , and M Arx7V are the sets of ./V-dimensional 
vectors of real numbers, of complex numbers, and of N x N 
Hermitian matrices, respectively; upper-case calligraphic let¬ 
ters and braces denote finite sets; E{-}, Cov[-], (-) T , (-)t, Z-, 
Tr [-], denote expectation, covariance, transpose, Hermitian, 
phase and matrix trace operators, respectively; Onxm (resp. 
In) denotes the N x M null (resp. identity) matrix; On 
(resp. ly) denotes the null (resp. ones) column vector of 
length N\ vec(M) stacks the first to the last column of 
the matrix M one under another to form a long vector; 
det(A) and ||A||i? denote the determinant and Frobenius 
norm of matrix A\ A® B indicates the Rronecker product 
between A and B matrices; diag(A, B) denotes the block- 
diagonal matrix obtained by placing matrices A and B along 
the main diagonal; the symbol means “distributed as”; 
x ~ denotes a complex (proper) Gaussian- 

distributed vector x with mean vector // £ C ;V x 1 and covari¬ 
ance matrix E £ C NxN ; X ~ CNnxm(A,B,C) denotes 
a complex (proper) Gaussian-distributed matrix X with mean 
A e C nxm and Cov[vec(X)] = B <g> C; S ~ CW N (K, A) 
denotes a complex central Wishart distributed matrix S with 
parameters K £ N and A £ C NxN positive definite matrix; 
M ~ CT a (A. £, rn) is a non-central multivariate complex F 
distributed matrix M with mean A and parameters a, £, and 
to; Pa denotes the orthogonal projection of the full-column- 
rank matrix A, that is Pa — [A(AtA) _1 At], while P\ its 
complement, that is P\ = (7 — Pa)- 

II. Problem Formulation 

We assume that a matrix of complex-valued samples X £ 
C NxK is collected, accounting for both primary (signal¬ 
bearing) and secondary (signal-free) data. The hypothesis 
testing problem under investigation can be formulated as: 

(Ho- X = A t B t C + No 

[Hi: X = (A t B t + A r B r ) C + N () ih 

where: 

• Nq £ C NxK is a matrix whose columns are independent 
and identically distributed (iid) proper complex normal 
random vectors with zero mean and (unknown) positive 
definite covariance matrix FI* £ C ,Vx,v , that is AT 0 ~ 

CNnxk(Onxk, Ik, R *); 

• B t £ C txM and B r £ C rxM denote the (unknown) 
deterministic matrix coordinates, representing the inter¬ 
ference and the useful signal, respectively; 

• A t £ C Nxt and A r £ C Nxr represent the (known) 
left subspace of the interference and the useful signal, 
respectively. The matrices A t and A,, are both assumed 
full-column rank, with their columns being linearly inde¬ 
pendent; 
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• Similarly, C £ C M x K is a known matrix describing the 
right subspace associated to both signal and interference; 
the matrix C is assumed full-row-rank. 


Additionally, aiming at a compact notation, we define A = 

_ n T 


''N x J 


and B = 


B B, 


£ C 


JxM 


, where 


we have denoted J = r + t. 

From inspection of Eq. |T|), we notice that the considered 
test has a complicated structure, which is thus difficult to ana¬ 
lyze. Therefore, before proceeding further, we first show that a 
simpler equivalent formulation of the considered problem can 
be obtained in the so-called canonical form ED- 

With this intent, we hist consider the QR-decomposition 
of A = Q a R a , where Q a £ C NxJ is a slice of a unitary 
matrix (i.e., Q),Q a = Ij) and R a £ C jxj a non-singular 
upper triangular matrix. It can be readily shown that Q a and 
Ra can be conveniently partitioned as: 


Qa 


\Qa,t 


Qa,r] 


Ra 


Ba,t Ba,x 

Orxi Ra : r 


( 2 ) 


where Q a t £ C Nxt and B at £ C txt arise from the QR- 
decomposition of At, namely A t = Q a ,t R a ,t > with Q a ,t 
such that Q lt t Q r ,.t = It and R,,.t a non-singular upper tri¬ 
angular matrix. Furthermore, £ C txr , and R r ,. r £ C rxr 
is another non-singular upper triangular matrix. Similarly, 
Qa.r £ C Nxr is such that Q^ a r Q a . r = It- Equalities in 
Eq. ([2ji are almost evident consequences of the well-known 
Gram-Schmidt procedure ED. Now, let us define a unitary 
matrix U a £ <C NxN whose first J columns are collectively 
equal to Q a . Then, it follows that: 




It ®txr 




A 4 

UlQa = 

0 rX £ Ir 

= [E t 

E r \ 

( 3 ) 


ec NxJ 

P(N-J)xt ®(N-J)xr_ 




where 

E t = 

[It 0 txr 0 tx(Ar _j)] T and 

E r 

A _ 


[0 rx t I r O rx (jv-j)] , respectively. Also, let C be ex¬ 
pressed in terms of its Singular Value Decomposition (SVD) 
as 

C = U 1 A 1 V?, (4) 


where U 1 £ C MxM and Vy £ C KxK are both unitary 
matrices, and the matrix of the singular values A 7 £ C MxK 
has the following noteworthy form: 


A 7 = 


■7 0 Mx(K-M) 


(5) 


with A 7 £ C MxA/ being a diagonal matrix. Therefore 

CV-y = M 7 [Im 0 m x (k-m)\ (6) 

holds, where Af 7 = C/ 7 A 7 . 

Given the aforementioned definitions, without loss of gen¬ 
erality we will consider the transformed data matrix Z = 
(UlXVy) £ C NxK in what follows. Such transformation 
does not alter the hypothesis testing problem being considered, 
as it simply applies left and right rotations to data matrix X 


(viz. multiplications by unitary matrices). The new data matrix, 
when Hi is in force, can be expressed as: 


Z = Ul (q q R a B t / 7 A 7 V 7 t) Vy + N 

= A ( R a B My) [. I M 0 mx(k-m)\ + N 


= A 


B t , i 

B 




J Mx 


(K-M)\ + N 


(7) 

( 8 ) 
(9) 


where we have defined B t i = (( R a ,t B t + R a ,x B r ) My) £ 
C txM , B = ( R ar B r My ) £ C rxM and N = 

(UlN 0 Vy) £ <C NxK , respectively. Furthermore, for the sake 
of notational convenience, we define B s = \B'j\ B 1 j . 
On the other hand, when Ho holds true, the matrix Z can be 
expressed as: 


Z = Ul [Q a R, 
B 


B t 

0 rxM 


Uy Ay Vj ) Vy + N 


= A 


•t, o 

OrxM 


[Im 0 Mx(K-M)\ + N 


( 10 ) 

( 11 ) 


where B t o = {R a .t B t My) £ C txM . Furthermore, aiming 
at keeping a compact notation, we will employ the definition 
C = [Im 0 M x(k-m)] i* 1 what follows. Gathering all the 
above results, the problem in Eq. 0 can be equivalently 
rewritten in terms of Z as: 


Ho: Z = A 


B t , o 
OrxM 


C + N 


Hu Z = AB S C + N 


( 12 ) 


Finally we recall that, since N 0 ~ CJ\f nxk( 0 nxk, Ik, R*), 
N is distributed as N ~ CM nxk(0nxk,Ik,R), where 

R±(uiR>u a ) m. 

An important remark is now in order. Specifically, for the 
problem in Eq. the relevant parameter to decide for the 
presence of a target is B,. Otherwise stated, if the hypothesis 
Hi holds true, then | B r \\ p > 0, while \\B r \\p = 0 under the 
target-absent hypothesis (Ho)- As a consequence, since R a , r 
is non-singular, problem in Eq. (fl2|) is equivalent to: 


Ho : ||B||f = 0 , 
Hi : 11.Blip > 0 , 


(13) 


which partitions the relevant-signal parameter space, say 0 r , 
as: 


0r = {OrxM}U {B € C rxM : 1|B|| f >0} . (14) 

©r,0 © r ,l 

The canonical form in Eq. ( [T2 ) i will be exploited hereinafter 
in our analysis. 

In the following, our analysis is carried out assuming that 
( K—M ) > N. Such condition is typically satisfied in practical 
adaptive detection setups (2). 


III. Maximal Invariant Statistic 

In what follows, we will search for decision rules sharing 
invariance with respect to those parameters (namely the nui¬ 
sance parameters, R, B t \, and B t ,o) which are irrelevant for 
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the specific decision problem. To this end, we resort to the so- 
called “Principle of Invariance’ ’ina, whose main idea consists 
in finding transformations that properly cluster data without 
altering 

• the formal structure of the hypothesis testing problem 
given by ( p~4| >; 

• the Gaussian assumption for the received data matrix 
under each hypothesis; 

• the double-subspace structure containing the useful signal 
components. 

The following subsection is thus devoted to the definition of 
a suitable group which fullfils the above requirements. 

A. Desired invariance properties 

Let 


• V(G, F) £ £, there exists (G_i,F_i) £ C such that 
(G_ 1 ,F_ 1 )o(G,F) = (G,F)o(G_ 1 ,F_ 1 ) = (G/,J» 
(Existence of Inverse element). 

Also, the aforementioned group leaves the hypothesis testing 
problem in Eq. ( fl2| invariant under the action I defined by: 

£{Z C ,S C ) = (GZ C + F,GS C G^ V(G,F)££. (22) 

The proof of the aforementioned statement is given in Ap¬ 
pendix [A] Moreover, it is important to point out that £ 
preserves the family of distributions, and, at the same time, 
includes those transformations which are relevant from a 
practical point of view, as they allow claiming the CFAR 
property (with respect to R and B t ,f) as a consequence of 
the invariance. 


V c ,i4 


J-M 
{K—M) x M _ 


v c , 2 ± 


*M x (K—M) 

Ik-m 


(15) 


Z c 4 (ZV c ,f) £ C NxM 
Z c ,± = (Z V c>2 ) £ C Nx{K ~ M) 


s c 4 (Z Ci± Zl ± ) = (ZP^tZt) £ c 


NxN 



r 

G n 

g 12 

g 13 ‘ 

5 = 

G = 

Orx£ 

g 22 

g 23 


l 

P(N-J)xt 

0(JV-J)xr 

g 33 


B. Derivation of the MIS 


In Sec. Ill-A we have identified a group £ which leaves 


and observe that P c t = (VcpV^f) and PA = (V Ci2 V^ 2 ). 

Also, let us consider the sufficient statistical { Z c , S c }, where 
the mentioned quantities are defined as 


(16) 

(17) 

(18) 


Clearly, given the simplified structure of G, Z c (resp. Z c ±) is 
simply obtained by taking the first M (resp. the last K — M) 
columns of the transformed data matrix Z. 

Now, denote by Q£(N) the linear group of N 
singular matrices and introduce the following sets 


unaltered the problem under investigation. It is thus reasonable 
finding decision rules that are invariant under £. Toward this 
goal, the Principle of Invariance is invoked because it allows 
to construct statistics that organize data into distinguishable 
equivalence classes. Such functions of the data are called 
Maximal Invariant Statistics and, given the group of trans¬ 
formations, every invariant test may be written as a function 
of the maximal invariant CD- 

Before presenting the explicit expression of the MIS, we 
give the following preliminary definitions based on the parti¬ 
tioning of matrices Z c and S, : 


(23) 


x N non- 

Zf 


'S n 

Sv 2 

Sis 

z c = 

Z 2 

; S c = 

S 2 i 

S 22 

S 23 


Z 3 


S 31 

S 3 2 

S 33 




£ Q£(N)( 19) 
G : Gn £ GC(t), G 22 £ Q£(r), G 33 £ S£(N - J)} 

£ C NxM : Pi £ C txM } (20) 


F = 


Pi 

OrxM 


(N-J)xM. 


along with the composition operator “o”, defined as: 


where Z x £ C txM , Z 2 £ C rxM , and Z 3 £ C 
respectively; S l3 , ( i,j ) £ {1,2,3} x {1,2,3}, is a sub-matrix 
whose dimensions can be obtained replacing 1, 2, and 3 with t, 
r, and (N — J ), respective[\0 Furthermore, we also define the 
following pardoning for Z Cj _l, which will be used throughout 
the manuscript: 


= [z 


T 

_L,1 


yT 


ZI >3 ] 


(24) 


(G a ,F a )o(G b7 F b ) = (G b G a ,G b F a + F b ) (21) 

The sets and the composition operator are here represented 
compactly as £ = {Q x P, o). Then, it is not difficult to 
show that £ constitutes a group, since it satisfies the following 
elementary axioms: 

• £ is closed with respect to the operation “o”, defined in 
Eq. {2T}; 

. V(G a , F a ), ( G b ,F b ), and (G C ,P C ) £ £: [(G a ,F a ) o 
(Gfc, F b )} o (Go, Pc) = (G q , F a ) O [(G b , F b ) O (Go, F c )\ 
(Associative property)', 

• there exists a unique (G/, P/) £ £ such that V(G, P) £ 
£■ (G u Fi) o (G,P) = (G,P) o (Gj,Fj) = (G,P) 
(Existence of Identity element ); 


where Z ± ,i £ C tx( - K ~ M \ Z L>2 £ C rx{ - K ~ M> > and Zj_, 3 £ 
q(n-j)x(k~m), respectively. Observe that each sub-matrix 
of S c in Eq. ( p3| can be expressed in terms of Eq. (24 1 , 
that is, Sij = (Z± t iZ]_ j). We are thus ready to present the 
proposition providing the expression of a maximal invariant 
for the problem at hand. 


Proposition 1. A MIS with respect to £ for the problem in 
Eq. ([72]) is given by: 


T^{Z1 3 SP 3 Z 2 ,} 


T(Z C , So) = <M T b 4 {z{ S'-, 1 Z 3 } 


Z\Sf 2 1 Z 2 


J < N 


J = N 


(25) 


'indeed, Fisher-Neyman factorization theorem ensures that deciding from 
{Z c , S c } is tantamount to deciding from raw data Z f3l . 


2 Hereinafter, in the case J = N, the “3-components” are no longer present 
in the partitioning. 
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Figure 1. Block diagram of processing leading to a generic invariant test. 


where Z 2 .3 = ( Z 2 - S 23 S 33 Z 3 ) and S 2.3 — (S 22 - 
S23 S 33 S32 ) ■ 

Proof: The proof is given in Appendix [B] ■ 

Some important remarks are now in order. 

• In the case J < N, the MIS is given by a pair of matrices 
(namely T a and Tf) where the second component (T/,) 
represents an ancillary part, that is, its distribution does 
not depend on the hypothesis in force; 

• In the case J < N, the matrices T a £ C MxM and Tf, £ 
C MxM have rank equal to min{M, r} and min{M, N — 
J}, respectively; 

• It is of certain interest comparing the general expression 
in Eq. (|25|) with the MIS instances obtained in lfl4l . 
m, mrm for specific adaptive detection scenarios. 
Accordingly, Sec. [V] will be devoted to comparisons and 
exhaustive discussion of the specialized forms in some 
relevant scenarios; 

• Finally, exploiting IflOl Thm. 6.2.1], every invariant test 
may be written as a function of Eq. ( |25j ) (see. Fig. |T] for a 
schematic representation). Therefore, it naturally follows 
that every CFAR test can be expressed in terms of the 
MIS. Part II of this study will be devoted to the design 
of theoretically-founded detectors whose CFARness will 
be proved by showing their dependence on the data solely 
through the obtained MIS. 


IV. Statistical Characterization of the MIS 


as in Eq. (231 and the pair ( G°,F° ) is suitably defined as: 


G° = 


F° = 


Gf 

G \2 

g ° 13 


Orxl 

7>— !/ 2 

■ rL 2.3 

n n-1 

; (28) 

_0 (N-J)xt 

0(A-J)xr 

D- 1 / 2 

r *'33 


F° 




0 rxM 

_0 (N-J)xM_ 

5 


(29) 


where G°j, i = 1, j £ {1,2,3}, and F° are generic matrices 
of proper dimensions, while R 22 £ C rxr , R 23 £ C rx ( N ~ J \ 
and R 33 £ C( N ~ J ' >X ( N - J ' > are obtained partitioning the 
true covariance matrix R in the same way as done for 
S in Eq. (23 1 . Finally, we have defined R 2.3 = R 22 — 
R23 R33 R' 23 - 

Hereinafter, we will study MIS statistical characterization 
after the trasformation (G°, F° ). This will simplify the sub¬ 
sequent analysis and does not affect the obtained results since 
the MIS is (by definition) invariant with respect to every 
trasformation belonging to C. 

Now observe that, under 7T;, i £ {0,1}, it holds: 


I Ut = G 3 


Z 2 

^3 


| 


' CJT (n- 


( N-t)xM 


i Go 


J 22 

*32 


J 23 

?° 

5 33 


= <G 


Ir 

0(AT- J)xr 
^22 


B, Im, J-N—t 


532 


S23 

S 33 


(Gl)' 


CW N ~t(K — M, I N _ t ) 


(30) 


(31) 


In this section, we provide the statistical characterization of 
the MIS for the case J < N and then, we will give a corollary 
referring to J = N. To this end, we show that the MIS can be 
written as a function of whitened random vectors and matrices 
and then we find a suitable stochastic representation by means 
of a one-to-one transformation. 

First, we consider the following transformation (G°, F°) £ 
C , which leads to: 


where we have debited G% £ C^ N t '> x ( N *) 


as 


G°3 = 


R 


-1/2 

2.3 


(N-J)xr 


TD p i 

-n 2 3 u 2 3 -K33 

>— 1/2 


R 


■33 


(32) 


Thus, exploiting the invariance property, we can equivalently 


rewrite the two components of T(Z C 1 S C ) (cf. Eq. (25 1 ) in 
terms of the whitened quantities 


T a = (Z 2 ° 3 )t (SI 


J 2.3’ 


(33) 


and 


Z° = G°Z C + F° = 


S° = G° S r G ot = 


zf 

Z°2 T 

ZfV 

(26) 

T b = (Z ° 3 )t (S^a)" 1 Z 3 °, 

(34) 






where Z 2 ° 3 4 - 5 2 ° 3 (^g)" 1 Z 3 ° and S° 2 . 3 

A QO 

— °22 

Sl-3 

S° 12 

Q° 1 
°13 



S 23 (S 33 )~ 1 S 32 , respectively. Let us focus on Z 23 

and rewrite 

C° 

*^21 

(f° 

°22 

°23 


(27) 

Z° 2 = [ Z 2,l ''' Z 2,m\ 

(35) 

*^31 

Q 0 

°32 

°33J 



(K-M) 



where Z° and S| (i,f, m £ {1,2,3}) are similarly debned 


s 2 ° 3 (s^r 1 z°3 = E r U is™r 1 z° 3 (36) 


k =1 
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^ 2.3 = [9i '' ■ Qm] = (* 2,1 ^ EkL i M r 2 ,k r l,k a i) ■ ■ ■ (* 2 ,m - Ek=i M r^ k r\ k a M ) 


(37) 


where r,,*, generically denotes the fc-th column of Z±j. Given 
the aforementioned definitions, we obtain the explicit form 
given in Eq. (37 1 (at the top of next page) for Zj 3 in terms of 
its columns qp, i = 1,... M, where we have further defined 
a,£ = | {S 33)- 1 z 3 and 2,3, £ similarly represents the £-th 
column of Z 3 , that is Z 3 = [z 3il • • • z 3 ,m] ■ 

First, we observe that r^,k ~ CAf r (O r , I r ) and r 3 ^ ~ 
CMn- j(On~j, In-j)', also it is apparent that these vectors 
are all mutually independent. Before proceeding, we define 
the short-hand notation “#3” to denote the conditioning with 
respect to all the terms with subscript “3”. Then, it can be 
shown that qi\(#3, Ho) is Gaussian distributed (recall that 
Z 2 i\Ho ~ CM r (O r , I r )) with mean vector 0 r and covariance: 


{ K-M K-M 

( z v ~ 51 r 2,fc r 3,fca/)(*2,«- 55 r 2* r i,k a *) ] 

k= 1 fe=1 

I r (l + (5 3 0 3 )- 1 23,f) (38) 

Similarly, the cross-covariance between 9«|(#3,Ho) and 
q m |(#3,7f 0 ) is given by: 



i K—M K-M \ 

( Z 2,m 55 r '^- r 5 a '")( Z V ^ 55 r 2 > feT'l ;fc af) t l = 

fc=l k= 1 J 

T ( z l, m (^3°3)- lz 3, f ) (39) 

Therefore, in view the aforementioned results, it follows that 
C 2.3 = vec(Zf 3 ) is conditionally distributed as: 

C2. 3 |(#3,fto) ~ CJVrM (0 rM , (Jm + (Z 3 °)t(5 3 ° 3 )- 1 Z 3 °) ® I r ) 

(40) 


components with subscript “3”). We then consider a one-to- 
one transformation of T(Z°, S°), defined as: 


Ti(Z° c , S ° c ) 4 


Kl (Z|. 3 )t (Sia )- 1 Z. 
(Z|)t (SS 3 )-i Z 3 ° 

(S 2 ° a)" 1 X 
(Z 3 °)t (S| 3 )-i Z| 


2.3 


K 


s3 


A_ 

Ti.a 


Ti ,b_ 


(44) 

(45) 


It is clear that, since Ti(-) is a one-to-one transformation of 
the MIS, it is a MIS itself fH. Therefore, without loss of gen¬ 
erality, we will concentrate on the statistical characterization 
°f Ti(-). 

We start by recalling that S% 3 is independent of { S % 2 , S 33 } 
l20l Thm. A. 11]. Also, we notice that S 33 ~ CY\Jn-j(K — 
M, I N _j) and SI 3 ~ CW r ((if - M) - {N - J), I r ). 
These results hold under both the hypotheses. Furthermore, 
conditioned on Ho, X is independent on 2j h (as it is 
independent on terms with subscript “3”). 

Therefore, it follows that conditioned on Ho, T-\ a and Jj h 
are statistically independent matrices, which means that the 
joint pdf can be written as fo(T 1>a ,T ltb ) = f 0 {T lta ) 

(as T-\ h denotes the ancillary part of the MIS and thus its pdf 
is independent on the specific hypothesis). Finally, it is worth 
noticing that in the case M < r we obtain the explicit pdf of 
T lt a\H 0 ~ CT m (0mxm, r, (K - M ) - (N - J)) and, for 
M <{N- J), T lih ~ CT m (0mxm, N-J, (. K - M) - 
(N — J )), following ETI . 

On the other hand, when Hi holds true, it is not difficult 
to show that Z% 3 \(Hi,#3) ~ CMrxM{Ro\ /2 B, (I M + 
(Z 3 0 )t(5 3 ° 3 )- 1 Z 3 0 ),/ r ) and consequently X\(H U #3) ~ 
CJVrxM (R-^BK^Im, I r ).A direct inspection of the 
last result reveals that 


Then, we whiten £2.3, that is, we define: 


A 

X — 


( I M + Z 3 0t (S 3 °; 


\-l 


337 


^3°) 


- 1/2 


^2.3 (41) 


CM 


rxM 


X\(#3,Hi) = X\(Ti yb ,Hi) ~ 

(r- 1 ^ 2 B ((J m + T, b )- 1/2 ) T , J M , J r ) 


(46) 


which evidently gives x\(#3,Ho) ~ CM r m (0 r m , IrM)■ 
Also, we observe that: 


(Jm + Z^iSl^zi) c 

(j M + z 3 ^(s , 33 )- 1 z 3 ) 


5 J 

” 1/2 


- 1/2 

5 Jr 


(42) 


which readily follows from the distributive property of Kro- 
necker product. As a consequence, we have that x = vec(X), 
where we have defined 


x = (Z |. 3 K s3 ) , 


(43) 


and K s3 4 


_(J M + Z 3 ot (S 3 ° 3 )- 

straightforward to show that X\(Ho,#3) = X\Hq ~ 
CM r x m (Orx m , Im, Ir) (i.e., it does not depend on the 


wzi )- 1 / 2 


Also, it is 


which underlines that Tj )Cl and T\j, are statistically de¬ 
pendent under Hi, thus leading to /i(Tj 0 ,Ti&) = 

/i(Tj a \Ti b) f(Tb). Again, in the specific case M < r it holds 
ri.„|(W 1 ,Ti, 6 ) ~ CT M (n, r, (.K - M) - (N - J)), where 
we have denoted $1 = (/T ( 3 Ji^g B K s3 ), following 12T1 . 

Finally, we conclude the section with a discussion on 
the induced maximal invariant in the parameter space ED- 
The induced maximal invariant represents the reduced set of 
unknown parameters on which the hypothesis testing in the 
invariant domain depends. It is not difficult to show that for I- 
GMANOVA model this equals to T p = B R.,l, B £ <C MxM . 
In addition, the induced maximal invariant is not full rank 
in the general case, with corresponding rank being equal to 
min{r, M}. It is worth remarking that such result applies in 
general, that is, the distribution of the MIS will depend on 
the parameter space only through T p , following classic results 
from HOl . 
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V. MIS IN SPECIAL CASES 


A. Adaptive detection of a point-like target 

In the present case we start from general formulation in 
Eq. Q and assume that: ( i) t = 0 (i.e., there is no interfer¬ 
ence); (ii) r = 1 (thus J = 1) i.e., the matrix A r collapses 
to a r £ C Nxl ; (in) M = 1, i.e. the matrix B r collapses to 
a scalar b r £ C and (iv) c = [l 0 • • • 0] £ C lxK (i.e. a 
row vector). Such case has been dealt in 1(121 . Therefore, the 
hypothesis testing in canonical form is given by: 


•Ho : Z = N 

Hi : Z = ab c + N 


(47) 


where a = [l 0 • ■ ■ 0] T £ C Nxl and c = c. By looking 

at the general MIS statistic expression in Eq. ( |25) l, it is not 
difficult to show that the present problem admits the following 
simplified pardoning: 


Therefore, the hypothesis testing in canonical form is given 
by: 


Ho : Z = N 
Hi: Z = Abc + N 


(55) 


where A = [l r 0, rx (jv- r )] T £ (^Nxr anc j c = c. By 
looking at the general MIS statistic expression in Eq. © 
it is not difficult to show that the present problem admits the 
following simplified pardoning: 


Z c ,± — 


Z± : 2 

Z ±,3 




S22 

S 3 2 


S 23 

S33 


(56) 


where z 2 £ C Jxl , z 3 £ Z ±:2 £ C Jx ^ K ~^, 

z x , 3 £ C ( JV - 7 ) X ^- 1 ), S22 £ C jxj , S 23 £ C JX ( N ~ J K 

S32 £ C ( - N ~ j ' )xj and S33 £ respectively. 

Exploiting the above partitioning, gives S2.3 £ C jxj and 
the simplified expression: 


z 


C 


Z2 

Z 3 


Z C ,_L = 


z±, 2 

3 


S c 


S 22 
^32 


8 23 

S 3 3 


, (48) 


where z 2 £ C, z 3 £ C( JV_1 ^ xl , z ± :2 £ C lx ^ K ~ 1 ' 1 (i.e., 
a row vector), Z ^.3 £ S 22 £ C, S 23 £ 

(Qix(tv— 1 ) (-j e a row vec tor), s 32 £ C^ -1 ^ 1 and S 33 £ 
C( w —i)x( Ar_ 1 ), respectively. Exploiting the above pardoning, 
gives the simplified expressions: 

s 2.3 = (^22 ^23 ^33 ^ 32 ) (49) 

Z2.3 = (Z 2 — ■ s 23 S 33 z 3 ) (50) 


*2.3 = (z 2 - S 23 Si 3 1 z 3 ) £ C Jxl (57) 

Since z 23 and z 3 are both column vectors, it is not difficult 
to show that the two components of the MIS are both scalar 
valued and equal to: 

la = z 2 3 S 2 3 Z 2 . 3 , t b = Z 3 S 33 Z 3 , (58) 

which is the classic result obtained in m. 

Finally, the (scalar-valued) induced maximal invariant 
equals t p = b' Rf\b which is the result obtained in M. 
being equal to the SINR. 


which are both scalar valued. Therefore, it is not difficult to 
show that the two components of the MIS are both scalar 
valued and equal to: 


t n = 


\Z2 -S 23 S 3 3 Z 3 1 :2 
[ s 22 — S 23 S 33 S32] 


(51) 

(52) 



Z2 — Z±pZ\^ 3 (Z l >3 Z\, 3 ) 1 z 3 

2 

Z_L,2 

Ik-i ~ Z1 3 (Z ± , 3 Z1 3 )-1 Z± , 3 

z ±,2 


t b = z\ H 33 1 z 3 

which can be rewritten in the more familiar form: 


t n . = 


t b = z f 3 (Z ±j3 Z [^) x z 3 

which are those obtained in nm 


(scalar-valued) induced maximal invariant is t p = (1| 2 / 7 * 2 . 3 ), 
where r 2 .3 = (r 2 2 — t'23 R 33 t'32), where we exploited 
similar simplified partitioning for R as for S c . The induced 
maximal invariant clearly coincides with the Signal-to-Noise 
plus Interference Ratio (SINR). 

B. Adaptive vector subspace detection 

In the present case we start from the general formulation in 
Eq .Q and assume that: (i)t = 0 (i.e. there is no interference, 
thus J = r); (ii) M = 1, i.e. the matrix B r collapses to a 
vector b r £ C Jxl and (Hi) c = [ 1 0 ••• 0] £ C lxK 
(i.e., a row vector). Such case has been dealt in nm M- 


C. Adaptive vector subspace detection with 
structured interference 

In the present case we start from general formulation in 
Eq. 0 and assume that: (i) Mj= 1, i.e., the matrices B r 
and B t collapse to the vectors b r £ C rxl and b t £ C txl , 
respectively; (ii) c = [l 0 0] £ C lxK (i.e. a 

row vector). Such case has been dealt in m. Given the 
aforementioned assumptions, the problem in canonical form 


(53) 

is given as: 

1 

Ho : 

Z = A 

(54) 

1 

Hi : 

1 _ 1 

II 

tsi 

the 

where b t ,i £ C txl , 

b £ C rxl 




ol 

b T 


c + N 

t (59) 

c+N 


and c = c. By looking at the 
general MIS statistic expression in Eq. ( |25j ) it is not difficult to 
show that the present problem admits the following simplified 
expression: 

z l = [zf z\ zj ] 1 (60) 

where zi £ C txl , z 2 £ C rxl and z 3 £ C ( N ~ J ) xl , 

respectively. Therefore, it is readily shown that S 2 . 3 £ C rxr 
and z 2 .3 = (z 2 - S 23 Sf 3 z 3 ) £ C rxl . Since z 2 .3 and z 3 
are both column vectors, it is not difficult to show that both 
components of the MIS are scalar-valued (similarly to the “no¬ 
interference” case) and equal to: 


tn = Z. 


1 Sr, * Z2.3 ; t b = zt s. 1 


2.3 ^2.3 


^33 


*3 5 


(61) 
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which can be recognized as the result obtained in IT8l . It is 
worth noticing that this result seems identical to that obtained 
in the previous sub-section (i.e., the interference-free case). 
However, we observe that, as opposed to the expression in 
definition of constituents of MIS in Eq. ( pH) is obtained 
by discarding terms with subscript “1”. In other terms, Eq. ( |6T) 
is analogous to ( |58) only after projection in the complementary 
subspace of the interference. 

Finally, the (scalar-valued) induced maximal invariant is 
t p = h' R 2 \ b, which coincides with the result obtained 
in Da, being equal to the SINR in the complementary 
interference subspace. 


D. Multidimensional signals 


In the present case we start from general formulation in 
Eq- 0 and assume that: ( i) t = 0 (i.e., there is no interference, 
thus J = r), (ii) A r = I jy (thus J = r = N) and (Hi) 
C = [Im Oat x(k-m)]- Such case has been dealt in lfT6l . 
Therefore, the hypothesis testing problem in canonical fornj^] 
is given by: 


Ho : Z = N 
Hu Z=BC+N 


(62) 


By looking at the general MIS statistic expression in Eq. ( |25) 
it is not difficult to show that the present problem admits the 
following simplified partitioning: 


z c — Z-2 , Z c ,_ I — Z±, 2 , Sc — S-22 , (63) 

where Z 2 £ C NxM , Z ^ 2 G C Nx ^ K ~ M) and S 22 G C NxN , 

respectively. Since in this particular setting J = N holds, 
we exploit the alternative expression for the MIS in Eq. ( [25) , 
which shows that the MIS reduces to a single matrix, being 
equal to: 

T{Z2,S22)=Z\S 22 1 Z2 

= z\ (Z^zlz )- 1 Z 2 (64) 


In the latter case, it is not difficult to show that the maximal 
invariant induced in the parameter space reduces to T p = 

B^R'B. 

It is now interesting to compare the result in Eq. (64 1 with 
that obtained in DU- Indeed, in the aforementioned work, the 
elementary action £(■) is defined as: 


i{Z 2 , S 22 ) = (g 22 z 2 U d , G22 S22 G 


(65) 


VG 22 e gc(N),\/U d eU{M). 

which, compared to Eq. ( [22) , enforces an additional invari¬ 
ance with respect to a right subspace rotation, via the unitary 
matrix U d ■ Clearly, this restricts further the class of invariant 
tests. Moreover, in Eq. (65 1 we have used U(M ) to denote the 
group of unitary M x M matrices. It was shown in fl6l that 
the MIS for the elementary action defined in Eq. © is given 
by the non-zero eigenvalues of the matrix 

- 1 / 2 (y 

2 V Z 2 


T C =S~ 


= (Z 


_L,2-< 


A 

J 3.,2> 


y\\ o -1 / 2 
* 22 

X 1/2 


(Z 2 Zl) (Z Xi2 Zl 2 ) 


t t-1/2 


( 66 ) 


denoted with eig(T c ) in what follows. Remarkably, we show 
hereinafter that the MIS in Eq. ( |66) can be directly linked to 
the expression in Eq. ( [64) . We first notice that, after defining 

Z m 4 |(Zj_ i2 Z} 2 )- 1 /2 Z 2 J £ C NxM , the following 
equalities hold: 

T = (Zl Z m ) T c = (Z m Zl) (67) 

Therefore, by construction, the matrices T € C MxM and 
T c £ C NxN are such that eig(T c ) = eig(T) holds (and the 
vector length equals min{M, N}), where we have expressed 
the non-zero eigenvalues through the implicit vector-valued 
function eig(-). Then, we notice that the action f 2 (-) can be 
re-interpreted as the composition of the following sub-actions: 

C 2 ,a {z 2 , S 22 ) = (G 22 z 2 , G 22 S 22 G\ 2 ) vg 22 e gc(N) 
hA^S 22 ) = {Z 2 U dl S 2 2) VU d £U(M). (68) 


It is then recognized that f 2 , a (’) = £(•) for the case of 
multidimensional signals. Previously, we have shown that the 
MIS for the elementary action £ 2 ,a(') = £(■) is simply given 
by the matrix T in Eq. (64 1 . 

Additionally, we notice that, for each U d £ U(M), 

T(Z 2 ,S 22 ) = T(Z 2 ,S 22 )=» 

T{Z 2 U d , S 22 ) = T(Z 2 U d , S 22 ) (69) 

Now, define the action £ 2b {-) as: 

l* 2 , b (T) = ( U\TU d ) \/U d £ U{M). (70) 

where T £ H MxM . It is not difficult to show that a MIS for 
the elementary operation b (-) in Eq. (70 1 is given by eig(T). 
Therefore, exploiting flOl p. 217, Thin. 6.2.2], it follows 
that a MIS for the action £ 2 (-) is the composite function 
eig(T(Z 2 , S 22 )). However, since as underlined in Eq. (67 1 , 
we have eig(T) = eig(T c ), this clearly coincides with the 
result in Da. 

Finally, by similar reasoning it is not difficult to show that, 
in such a case, the induced maximal invariant is given by 
eig(Tp) = eig (B^R^B) = eig(J?- 1 / 2 B B^R- 1 / 2 ), thus 
obtaining the result in Da. 


E. Range-spread Targets 

In the present case we start from general formulation in 
Eq. (|T) and assume that: (i) t = 0 (i.e., there is no interference, 
thus J = r); (ii) r = 1 , thus the matrices A r and B r 
collapse to aT r £ C^ xl and b r £ C lxM (i.e. a row vector), 
respectively; (Hi) C = [Im OmxA-m]- Such case has 
been dealt in E2, D3- Therefore, the hypothesis testing in 
canonical form is given by: 


H 0 : Z = N 
Hi_: Z = abC + N 


(71) 


where a = [l 0 ■ • • 0] T £ C Nx \ b £ C lxM and G = 

G, respectively. By looking at the general MIS expression in 
Eq. ([25}, it is not difficult to show that the present problem 
admits the following simplified partitioning: 


3 It is worth noticing that in this case the original formulation in Eq. [Tj is 
in canonical form already. 




-Z2 

Z 3 


Z c ,± — 


Z -L,2 

Z±,3 


S r = 


S 22 ^23 
•S 32 S 33 


(72) 
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where 2 2 G C lxM (i.e., a row vector), Z 3 G C( A,_1 ) xM , 
Zj _ i2 G C(i.e., a row vector), Zj _ j3 G 

(j^(w-i)x(A-M)^ g22 g C, g23 g (Af— 1 ) (j e ^ a row 

vector), s 32 G C ^ -1 ^ 1 and 533 G re¬ 

spectively. Exploiting the above partitioning, gives S 2.3 = 
(S 22 — •S 23 Sf 3 S 32 ) G C (i.e., a scalar) and the simplified 
expression: 

22.3 = («2 - «23 S 3 i Z 3 ) G C lxM (73) 

Given the simplified expressions for z 2 .3 (row vector) and 
S 2.3 (scalar), it is not difficult to show that the two matrix 
components of the MIS are given by: 

Ta = (—) 4 3 * 2.3 T b = z\ S 3 - 3 1 Z 3 (74) 

V S 2.3 / 

where the matrix T a is rank-one in this specific case (as it is 
the output of a dyadic product). Also, the induced maximal 
invariant in the parameter space equals T p = (^T_) b b. i.e., 
a rank-one matrix. 

It is now of interest comparing the MIS represented by 
Eq. ( f74) > with that obtained in fl7l . The approach taken in the 
following is similar to that used for multidimensional signals 
in Sec. |V-D| However, due to the more tedious mathematics 
involved, we confine the proof to Appendix [C] and we only 
state the results hereinafter. 

Indeed, in the aforementioned work, the elementary action 
£(•) is defined as: 


V Ud G U(M), where T b G H MxM and T a = (aal) (that is, 
a rank-one matrix). It is shown in Appendix JCl that the MIS 
for the elementary operation £ 2b {-) ' n Eq. (791 is given by 
{eig(Tb), eig(T a + T b )}. Therefore, by exploiting (TOJ P- 217, 
Thm. 6.2.2], it follows that the MIS for the action £ 2 (-) is the 
composite function 


UlzlSiZt) 

\eig S33 1 Z 3 + (^-) 4 3 *2. 3 ) 

which clearly coincides with the result in fl7l . 

Finally, it is not difficult to show that the induced maximal 
invariant in such a case can be obtained as eig(T p ) = (yyy) = 
||6|| 2 (at Il ] a) (since the rank-one induced maximal invari¬ 
ant has only one non-zero eigenvalue), which represents the 
overall SINR over the M cells, as defined in fTTM . 


F. Standard GMANOVA 

Finally, in the present case we start from general formulation 
in Eq. O and assume that: (i) t = 0 (i.e., there is no 
interference, thus J = r). Such model clearly coincides with 
that analyzed in 12, 0, unfortunately not dealing with the 
derivation of the MIS. Therefore, the hypothesis testing in 
canonical form is given by: 

h-- z=N (8D 

\Hi: Z = A B C + N 


e 2 (Z c , s c ) = (G Z c u d , G S c Gt) , (75) 

VGGg, VU d £ U(M), 

which, compared to Eq. ( [22] ), enforces an additional invariance 
with respect to a right subspace rotation of primary data, via 
the unitary matrix Ud- It was shown in lfl7ll that the MIS 
for the elementary action defined in Eq. ( [75] ) is given by the 
eigenvalues of the matrices 


(T a +T b ), T b , (76) 


denoted with eig(T a + T b ) and eig(T b ) in what follows. 
Remarkably, we show hereinafter that the MIS in Eq. (76 1 can 
be directly linked to the expression in Eq. We first notice 
that the action l 2 {-) can be re-interpreted as the composition 
of the following sub-actions: 


h,a{Z c , &c) = (G Z c , G S c G t ), VGGg, 
i2,b(Z c ,S c ) = (Z c U d ,S c ), VU d £ U{M). (77) 


It is then recognized that £ 2 ,a{') = £{') for the case of range- 
spread signals. Also, we have previously shown that a MIS 
for the elementary action £ 2 , 0 .(') = £(■) is given by Eq. (74 1 . 
Additionally, we notice that, for each U d G U(M), 


T{Z C ,S C ) = T(Z C1 S C )=> 

T(Z C U d , S c ) = T(Z C U d , S e ) (78) 


Now, define the action b (-) as: 

£* 2 , b {T a ,T b )=(u\ T a U d , U\ T b U d ) , (79) 


where A = [ij 0 ,jx(n-j)] T an d B G C JxM , respectively. 
By looking at the general MIS statistic expression in Eq. ( |25j ) 
it is not difficult to show that the present problem admits the 
following simplified partitioning: 


Z r .= 


Z 2 

Z 3 


Z r 1 = 


J J -,2 

?X,3. 


Sc = 


^22 

5 32 


^23 

^33 


, (82) 


where Z 2 G C JxM , Z 3 G Z ±j2 G <C Jx ^ K ~ M \ 

Z_ l, 3 G C ( n -J)x(K-M), g 22 g C JxJ S ’ 23 g C Jx(JV-J) ; 

S 32 G 0 N ~ J ' >XJ and 533 G C^ N ~ J ^ X ^ N ~ J \ respectively. 
Given the simplified definitions in Eq. ( [82] ), the MIS is readily 
obtained via the standard formula in Eq. <[25]). 

Finally, the induced maximal invariant is obtained through 
the standard formula T p = (£> '' R 2 ! 13 ). The sole difference 
consists in the rank of matrix T p , being equal to min{ J, M}, 
i.e., there is no reduction in the observation space due to 
structured interference. 


VI. Conclusions 

In the first part of this work, we have studied a generaliza¬ 
tion of GMANOVA model (denoted as I-GMANOVA) which 
comprises additional (deterministic) structured interference, 
modeling possible jamming interference. The study has been 
conducted with the help of the statistical theory of invariance. 
For the present problem, the group of trasformations leaving 
the hypothesis testing problem invariant was derived, thus 
allowing identification of trasformations which enforce CFAR- 
ness. Then, a MIS was derived for the aforementioned group, 
thus explicitly underlining the basic structure of a generic 
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CFAR receiver (several examples of CFAR receivers, based 
on theoretically-founded criteria, will be derived in part II of 
the present work). 

Furthermore, a statistical characterization of the considered 
MIS under both hypotheses was obtained, thus allowing for an 
efficient stochastic representation. As a byproduct, the general 
form of the induced maximal invariant in the parameter space 
was obtained for the considered hypothesis testing. Finally, 
the general MIS expression was particularized and compared 
with MIS obtained in specific instances found in the open 
literature. Analogies to other expressions of the MIS, obtained 
by enforcing invariance to a wider class of transformations (cf. 
Sec. |V-E| and |V-D| >, were underlined and discussed. 


Appendix A 

Invariance of the problem 

WITH RESPECT TO THE GROUP C 


In this appendix we prove the invariance of the hypothesis 
testing problem in ( |~i~ 2 j ) with respect to the group of trasfor- 
mations C defined in Sec. III-A Let ( G , F) £ C and observe 
that, tinder 77, i, the columns of GZ, + F are independent 
complex normal vectors with covariance matrix GRG ^ and 
mean: 


GAB S + F = 


GnB t i + G 12 B + F\ 

G 22 B 


0 


(N-J)xM 


B t , 1 

B 

.0 (N-J)xM_ 


= ab q 


(83) 

(84) 


where we have employed the definitions B tl = (< G\\B t 4 + 
G 12 B + i*ii) € C txM and B 4 ( G 22 B) £ C rxM , 
respectively. Also, aiming at compact notation, we have de¬ 
noted B s = [(B' t (B' ) t ] . Furthermore, it is not 

difficult to show that GS C G t = (GZ c> ±)(GZ Cj ±)^, with 
(G 2 C] i) ~ A/c( 0 jv x (a:-m)! Ik-m,G RG^). 

On the other hand, when 77o holds true, GZ r + F shares 
the same covariance structure as in the case of Hi, except for 
the mean, which becomes 



Bt, 0 
OrxM 


GuBt'0 + Fi 

GA 

+ F = 

0 rxM 

0 (N-J)xM 






OrxM 


(85) 

( 86 ) 


where B' t0 4 (GnB tfi + FJ £ C txM . Again, 
it is not difficult to show that (GZ C , j_) ~ 

CJ\f n x(k-m){0nx(k-m), Ik-m , GRG t). 

Therefore, it is apparent that the original partition of the 
parameter space, the data distribution, and the stmcture of 
the subspace containing the useful signal components are pre¬ 
served after the transformation ( G,F). Indeed, the following 
equivalence holds between the original and the trasformed test: 


Ho : ||B||jr=0 «=» ||B'|| f = 0, 

Hi : ||B|| f > 0 «=► ||b'|| f > 0, 


(87) 


where the nuisance parameters in the transformed space are 
B' t i and (GRG^). 


Appendix B 

Derivation of the Maximal Invariant Statistic 
In the present appendix we provide a proof for Prop. |T| 
In particular, hereinafter we will focus on the case J < N, 
as to the derivation for J = N can be obtained through 
identical steps. Before proceeding further, we recall that a 
statistic T(Z C ,S C ) is said to be a maximal invariant with 
respect to the group of transformations C iff 

(a) T(Z C ,S C ) =T[£(Z C ,S C )}, We£; ( 88 ) 

(b) T(Z C ,S C ) = T(Z C ,S C )=> 

3 l€C : (Z C ,S C )=£(Z C ,S C ). (89) 


Conditions (a) and (b) correspond to the so-called invariance 
and maximally properties, respectively. In order to prove (a), 
we first consider the following partitioning of matrix G and 
sub-matrix of S c : 


G = 


G 1 

0(AT-t)xt 


g 2 

Go ’ 




where S 2 £ COV-Oxfiv-t), d = 
[G 12 G 13 ] £ C ,x ^-‘) and 


S 2 2 

S 3 2 

Gn 


S 23 

$33 


(90) 


£ C txt , G 2 = 


Go 


A 


G 22 

0 (JV-J)xr 


G 20 

Goo 


g <£(tV-t) x (A-t) 


(91) 


Then, let (Z c , S c ) = £(Z C , S c ), with 


Z C = GZ C + F , S C = GS C G'. (92) 


It is apparent that the following equalities hold, when exploit¬ 
ing the specific structure of G and F: 

Z 2 = G 22 Z 2 + G 20 Zo, (93) 

Z 3 = Goo Zo , (94) 


and 


S 2 = Go S 2 G 


t 

3 > 


(95) 


where S 2 is similarly defined as in Eq. (90 1 . From Eq. (95 1 , 
it can be inferred that: 


’22 

— [G 22 G 23 ] 

S 2 [G 22 

G 23 ] * 

(96) 

’23 

= G 22 S 23 GL 

+ G 23 S 33 

<4 

(97) 

’33 

= G 33 Soo G 33 



(98) 


Additionally, exploiting the appropriate substitutions, it can be 
shown that: 


Z 2 .3 — (Z 2 — S 20 S 33 Z 3 ) — G -22 Z 2.3 (99) 

S2.0 = {S22 - S 23 S33 1 S02) = G22 S2.0 G \ 2 (100) 

Finally, substituting Eqs. ([98), and (QOOf into ([25), 
we obtain: 


T(£(Z C1 S C )) = 


’ y\ t 

^ 2.3 ^22 


(cy - 1 s^i g 22 


zt G \ 3 (G\ 3)- 1 53-3 1 (G 


r 

.3 ^ 2.3 
' Zo 

1 G 22 Z 

1 (r-t-1 


Z2 3 S2 3 Z2.3 

zl - 


33 


2.3 
) G 33 Z 3 


Z\ 3 ^2 3 ^2.3 

. ZlS^Zo . 


( 101 ) 

( 102 ) 

(103) 
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which thus proves (a). 

Now, in order to prove (b), assume that: 


T{Z C ,S C ) = T(Z C , S c 


" yt c-i 

't 


Z 2 3 S 2 3 Z 2 3 


^ S 3 3 ^3 


' c -1 

't o-i 


Z^2 3 S 2 3 -^ 2.3 


Z\S^Z 3 


(104) 

(105) 


The last equality can be recast as the following pair of 
equalities 


^3 4 = ^3^, 


>T 


y 3 y/ = Y 3 Y 3 t , 




(106) 


.-i/ 2 , 


J 2.3 


)t, Y 2 .3 4 (S ^ 1 / 2 Z 2 . 3 )\ Yj 4 


where Y 2 . 3 = (S 2 _ 3 

(S 3 ^ 2 Z^ and Yj = (S :i3 ! ' 2 Z 3 4■ It follows from direct 
inspection of Eq. (106i that there exist unitary matrice^] 

U 2 . 3 G C rxr and U _ 3 G C^ N ~ J '> x( - N - J '> such that Y 2 . 3 = 
Y 2 . 3 f / 2 . 3 and Y 3 = Y 3 t/ 3 . 

First, let us define the following block-triangular decompo- 

2 = L\ L 2 


sitions for matrices S 2 = L\ L 2 and S> = L\L 2 , where: 


Lo 4 


L 2 = 


4 

? l /2 

’2.3 


0 


o 1/2 o 
*33 *32 

el/2 

J 2.3 

q— 1/2 o 

*33 *32 


rx(N-J) 

? l /2 

5 33 


>rx(N-J ) 
°33 


(107) 


(108) 


Therefore, given the aforementioned definitions, it can be 
shown that: 


rt 

'2.3 


ryt 

1 2.3 

yt 

yt 

1 2.3 

t 






Zf\ 


(L\)~ l Z 2 3 = 

(109) 

u b (p 2 r x z 2 3 

; (HO) 

and U\ 4 

diag(t7+. 3 , C/t), 


where Z 23 

respectively. From comparison of right hand side, of Eqs. ( | 1 09[ ) 
and ( 1 10 [>, it readily follows that 


Z 23 = LlU 1 {L\)~ 1 Z, 


23 ■ 


( 111 ) 


From inspection of Eq. (jTTTJ, it is apparent that selecting 


the transformation G 3 = L ' 2 U\ (X/ 2 )" 1 (which is block 
triangular as dictated by Eq. ( |9T] i) automatically verifies the 
set of equations: 


G i Zyi = Z 


23 


g 3 s 2 g\ = s 2 


( 112 ) 


since it also holds 




L\)U^L 


l )- 1 


L\L 2 {L 2 ) 


2 

-irrt 


U{L 2 = 


U[L 2 = v 2 l 2 = s 2 


(113) 


4 Such property can be verified as follows: given the equality AA^ = BB^ 
between two generic matrices A and B and, after defining the eigenvalue 
decompositions (AA^) = UaAaU^ and (BB^) = UbAbU 
and the SVDs A = Ua X )a Vj and B = Ub X)# Vg, it is appar¬ 
ent that such equality implies: (i) Y>a = S 5 D; (ii) Ua = Ub D, 
where D and D denote diagonal matrices of phasors (recall that = 
X^XI^and Ab = SgSp. Thus, it follows from substitution that A = 
UbD Y>b D = Ub^b D D vj. Therefore, after defining the unitary 
matrix U* = Va (-D D)*Ug, it is finally demonstrated that AU* = B. 


Finally, we shown how to build the remaining blocks of G. To 
this end, let us denote S 3 = [S '12 S 13 ], S 3 = [S 12 S 13 ], 

Si = S11, and Si = Sn, and consider the block-triangular 
decompositions for matrices S c = L\ L c and S c = L\ L c as: 


L c 

L c 


A 


A 



s 1 / 2 


(Lir'sl 


ct 

*3 

0 


OfxGV-t) 

(u\l 2 ) 

tx(N-t) 

l 2 


(114) 

(115) 


where S.i = Si - S 3 S 2 " 1 S^ G C txt and S.i 4 
Si — S 3 S 2 1 S 3 G C txt . Also, since we need to ensure 
(G S C G 1) = S c , it suffices that 


L c Gl = L c 


- 1 

.C^I 

to 

Ot.x(N-t) 

6^1 0 tX (jv-t) 

(Lir 1 * 

% l 2 

G\ g\ _ 


sT 

0 tx(N-t) 



(Lilli )- 1 Si 

(U\l 2 ) 



(116) 


(117) 


from which the following set of independent equations arises: 


(0 sTg\ = sT 

(«) (.L\)- 1 SIg\+L 2 G\= (Llu^Sl 


which provides the “completing” solutions for matrix G : 

Gi = S 1 / 2 S 3 1/2 (119) 

G\ = S 2 X (g^sI - S^gI) (120) 

Up to now, we have shown how matrices Gi can be con¬ 
structed. Finally, it can be easily shown that matrix F-\ should 
be chosen as F\ = Z\ — ]U . ! =1 G 1 ,Y,. r I’his concludes proof 
for (b). 


Appendix C 

MIS OBTAINED BY ENFORCING 
ADDITIONAL INVARIANCE IN RANGE-SPREAD CASE 


In this appendix we show that the statistic 


xi = eig(T 6 ) 
x 2 = eig (T b + T a ) 


( 121 ) 


where T a = {aaJ) (a G C Mxl ) and T b G H MxM , is a MIS 
for the elementary action 

%, b (T a , T b ) = (U\ T a U d , U\ T b U d ) , (122) 

where U d G U(M). First, we observe that Eq. (121 1 is in 
one-to-one mapping with: 


(123) 


\xi= eig (T b ) 

\x 2 = \k\ 

where k = (U b a) and the modulus | 
be intended element-wise. Also, U b denotes the eigenvector 
matrix of T b , that is T b = U b K b Ul. The existence of the 
aforementioned mapping can be proved as follows. We start 


in Eq. (123 1 should 
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by observing that eig(T b + aaf) can be obtained as the zeros 
(with respect to the variable s) of the rational function ll23l : 

w(s) = (l + A:* (A b - s I M ) _1 fe) . (124) 

Also, since (A b — sIm ) _1 is a diagonal matrix, w(s) depends 
only on |fe|. Therefore eig(T), + aw) can be obtained starting 
from A b (viz. eig(T b )) and |fe|. Vice versa, the vector \x\ 
is obtained from eig(T b + an') and eig(T b ) by inverting 
Eq. ( |124[ ), that is: 


k) 

(A b — 

X 2 ,i 

Im) 1 k = —1, 


(125) 


M 

E 

n= 1 

(A b , 

\kn\ 2 _ 1 

n ^2 ,i ) 


(126) 




afe = -l, 

,M} 

(127) 

where 

^ b,n IS 

the 

n-th diagonal element of A b 

and 


€ 

A[|fci 

I 2 

••• \k M \ 2 ] T , 


(128) 

Oil 


,1 _ 

X 2 ,i)~ 1 ■■■ (>^b,N - X 2 ,i)' 

E T - 

(129) 


It is shown hereinafter that the linear system in Eq. (127 1 (with 
respect to the unknown vector e) admits a unique solution. 

Indeed, the generic ag represents a scaled version of 
{E p v a+bi ), where E p = diagje} and v a+b , denotes the i-th 
eigenvector of T b + a a |23, Eq. (2.1)]. However, since we 
assume that the eigenvalues are distinct with probability one, 
the eigenvectors v b j will be linearly independent. Therefore, 
it follows that also the set {oti}^L 1 constitutes a linearly 
independent basis. Such conclusion clearly implies that the 
system is invertible and admits a unique solution; therefore 
there exists a one-to-one correspondence between the statistics 
in Eq. ( |121 1 and ( 123[ >. 


Once established the correspondence between Eqs. (121) 


and (123), it suffices to show that Eq. ( | 123b is a MIS for 
the group of trasformations specified in Eq. ( | 1 22] ). In order to 
accomplish this task, we first prove invariance of statistic in 
Eq. ( |123[ ). Indeed, given the transformations: 

{f a = {U\aa)U d ), T b = (U\T b U d ), 


(130) 

It is readily shown that eig (T b ) can obtained as the zeros of: 

det (sJ - U ] d T b U d )=0 ^ det (si -T b )= 0 (131) 

thus coinciding with eig(T b ). Also, it holds 


\Ul'a\ = \UlU d U\a\ 
= \Ula\ 


(132) 


Therefore the statistic in Eq. (123) is invariant. We then prove 
maximality. Under the assumption 


eig (T b ) = eig (T b ) 
\U b a\ = \Ul a\ 


(133) 


it can be readily shown that there exists a unitary matrix 
V that ensures the equality {V''T b V) = T b , namely V = 
{U h D h Ul ), where D b is a diagonal matrix of arbitrary 
phasors. Similarly we have employed the eigendecomposi- 
tion T b = (Ui, A/, Ul). Additionally, in order to complete 


maximality proof, we need to prove that the aforementioned 
transformation, when applied to T a = an ', can be adjusted 
to satisfy: 


VfiaafiV = dal (134) 


After substitution, such condition can be rewritten as: 


U b D\ ul(aa^)U b D b u\ = aa* 

D\ Ul(aa))U b D b = ul&aW b 



D\ 


([/t°)] t =(t/ 6 t a) (f/^d) 1 


(135) 

(136) 

(137) 


The above rank-one matrix equality can be achieved by 
enforcing the vector equality 



(Ula)} =(ula) 


(138) 


by choosing each element of the diagonal matrix Dl in 
order to rotate each phase term of (Ul a) aiming at imposing 
Z (ul a^j = Z (Ula)j , since |[/ b a| = |£/ b a| by definition 
(cf. Eq. ( | 1 33| >). Therefore Eq. ( |123| ) (resp. Eq. ( | 1211 >) is a MIS 
for the aforementioned group of trasformations. 
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